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We report on a numerical experiment in which we use time-dependent potentials to braid non- 
abelian quasiparticles. We consider lattice bosons in a uniform magnetic field within the fractional 
quantum Hall regime, where v, the ratio of particles to flux quanta, is near 1/2, 1 or 3/2. We intro- 
duce time-dependent potentials which move quasiparticle excitations around one another, explicitly 
simulating a braiding operation which could implement part of a gate in a quantum computation. 
We find that different braids do not commute for v near 1 and 3/2, with Berry matrices respectively 
consistent with Ising and Fibonacci anyons. Near v = 1/2, the braids commute. 
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When two identical quantum mechanical particles ex- 
change places, the wavefunction typically acquires a 
phase: 9 = for bosons, and 9 = ir for fermions. Re- 
markably, there exist 2d systems whose "anyon" ex- 
citations display fractional statistics, with 9 0, 7T. Even 
more remarkably, there are models in which exchanging 
quasiparticles not only produces a phase, but also ro- 
tates the system between degenerate states [jl-Hi^ . Under 
these circumstances, exchanges may not commute. Ki- 
taev (l7j proposed using such nonabelian quasiparticles 
for quantum computation, with qubits constructed from 
the degenerate states. Quantum gates are implemented 
by "braiding" the quasiparticles: using timc-dcpcndcnt 
potentials to drag the quasiparticles around one another, 
switching their positions. The collective nature of the en- 
coded quantum information provides protection against 
various decoherence mechanisms. Here we start from a 
microscopic Hamiltonian, and numerically calculate the 
result of such a braiding experiment. We find that even 
for surprisingly small systems (4x4 lattices), this proce- 
dure can be used to establish non-abclian statistics, and 
hence to implement quantum gates. 

Explicitly calculating the results of a braiding opera- 
tion for a realistic microscopic Hamiltonian is difficult. 
Previous studies have focused on the properties of vari- 
ational wavefunctions in the limit of all quasiparticles 
asymptotically far apart (separation large compared to 
the magnetic length) [f| 0, EH-HH • As in physical exper- 
iments, a numerical experiment must contend with finite 
size effects, mixing of higher bands, the location of un- 
pinned quasiparticles, and uncertainty about both the ex- 
act many-body wavefunction and the interaction between 
a quasiparticle and the applied perturbation. Overcom- 
ing these difficulties is well worth the effort, since observ- 
ing the braiding of two quasiparticles provides a defini- 
tive test of exchange statistics. This numerical approach 
complements more indirect methods in real experiments, 
such as observing shot noise or interference effects in the 
tunneling of edge states [24 1 . 



We choose a model which is both experimentally rele- 
vant, and computationally tractable: hard-core bosons 
hopping on a square lattice, with phases on the hop- 
ping matrix elements corresponding to a uniform mag- 
netic field. This model describes Cooper pairs hop ping 
on a Josephson junction array in a magnetic field [25l — 
27| when the charging energy is large compared to the 
hopping energy. It also describes cold atoms in a deep 
optical lattice [28| with an artificial gauge field 291 33|. 
Recent developments in cold atom physics [33| suggest 
that the fractional quantum Hall regime will be attained 
in the near future. 

A general Hamiltonian for lattice bosons is 

// X [ J 3kz lhk a]a k + H.C.) ■ '-f-Y. a ] a ] a i a i W 
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a\/ctk creates/annihilates a boson at complex coordi- 
nate Zk on a square lattice with unit lattice spacing. 
Defining z = Zj — z^ = x + iy as a complex inte- 
ger, icf)jk = — ^ ( z j z * ~ z j z ) is the Pcicrls phase of 
the B field (with <f> the density of flux quanta per pla- 
quette). The properties of this Hamiltonian depend on 
the form of Jy. The simplest model would just in- 
clude nearest neighbor hopping (HI. As argued in [35| . 
the fractional quantum Hall states are particularly ro- 
bust if we use a particular gaussian hopping, Jjk = 

J(z) = J G(z)exp (-§ (1 - <)>) |z| 2 ) , where G (z) = 
(— ij + x +y+ x y anc | j Q i s a constant. We mostly take the 
hard-core limit of U2 — > 00. We define Jnn = Jo e_7r ^ 4 a s 
the energy scale of the problem. U3 is an artificial three- 
body repulsion which we introduce in some calculations. 
We showed in [35| that the single-particle spectrum of 
(Q]) reproduces the continuum lowest Landau level (LLL) 
with <fiL 2 degenerate single particle ground states on an 
L x L lattice. As explained in [35|, the longer range hop- 
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pings can be engineered by appropriately shunting the 
Joscphson junction array, or by appropriately tailoring 
the optical lattice potential. For <j> < 1/3 it suffices to 
include next-nearest-neighbor hopping. Since the lowest 
Landau level is preserved in ([T]), a LLL-projected calcu- 
lation in the continuum whould give similar results, at 
the cost of more complexity in the calculation. 

We add to Eq. (TTJ a time-dependent potential Vj(t), 
corresponding to a Hamiltonian H p = Y^j Vj{t) a ] a j- At 
time t = 0, we take V to be zero except on a few sites, 
where it is positive. We slowly change V such that 
Vj(T) = Vj(0), but with two of the potential bumps 
exchanged. If quasiparticlcs are pinned to the defects, 
this will exchange them. Experimentally, the potential 
Vj could be engineered by gates on individual Joscphson 
junctions, or through targeted lasers in an optical lat- 
tice. Such addressability was recently demonstrated in 
j36j . In our numerics we move our bumps by linearly 
reducing the amplitude of V on one site, while linearly 
increasing it on a neighbor. 

Under an adiabatic cyclic change of the Hamiltonian, 
non-degenerate states will return to themselves with an 
additional phase factor, while degenerate states can mix: 
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i f Edt Y, j M ij \ip j }. Throughout we ne- 



glect the J E dt term, where E{t) is the instantaneous 
energy at time t. This temporal phase can be experimen- 
tally distinguished from the geometric phase by travers- 
ing the path at different rates. The unitary matrix My 
is calculated by integrating the Berry connection: 



(2) 



M = Pexp 2?ri <p dA7 



Here, = i {ipi\ Va \ipj) is the Berry connection matrix, 
the \ipi) are a basis of degenerate states, A parametrizes 
the path, and P is the path ordering symbol. While the 
Berry connection 7 is a gauge-dependent quantity, the 
matrix M is physical and gauge invariant (up to joint 
choice of basis at the start and end points) . 

To numerically calculate Eq. ([2]), we use a method de- 



scribed in }37l |38|, breaking the path into many small 
discrete steps, engineered to maintain the degeneracies 
of the spectrum. For each point A on the path, we diag- 
onalize H to produce a basis ^(A)). This basis is not 
unique: the phases of |?/>i(A)) are arbitrary, but one can 
form a new basis by taking arbitrary linear superposi- 
tions of degenerate states. We fix this arbitrariness by 
choosing {ipi (A) \tpj (A + dA)) = 
matrix is then 



-O (dA 2 ) . The Berry 
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Following [38(, we generate the states \if>i(X + dA)) = 
\tpi(X)} = J2j{A~ 1 )ij\ipj(\)) by first determining the 

eigenstates |V>i(A + dA)) using a generic diagonaliza- 
tion algorithm, and then calculating the overlap matrix 



{tpi (A) \ipj (A+dA)) . Since A will be unitary only up 



to corrections of order dA, we perform a Gram-Schmidt 
orthogonalization at each step. 

In Fig. [I] we illustrate the initial configurations of the 
impurities and some of paths over which we move them. 
We use relatively small systems: between 3 and 9 par- 
ticles on lattices of up to 24 sites with periodic bound- 
ary conditions; with the hard core constraint the largest 
Hilbcrt spaces studied contained about 50,000 states. 
While state-of-the art algorithms on high performance 
computers would allow us to study larger systems, we 
find that finite size effects are already sufficiently small 
on these modest grids, presumably due to the robust na- 
ture of the topological effects of interest. Our algorithm 
was implemented in Mathematica on a desktop computer. 



RESULTS 

The results of our braiding calculations are summa- 
rized in table I. In all cases, the applied impurity poten- 
tials are strong. We assign each state an effective filling 
fraction i/ D g — -/V/-/Vlll, where A^lll is the number of 
single particle states in the LLL in the presence of the 
impurities. In every case studied, for N[ mp impurities 
TVlll = N$ — N[ mp (where is the number of flux 
quanta), showing that a full quasihole (QH) is pinned at 
each impurity, and in the thermodynamic limit v c f[ — > v. 
Each quasihole is a first order zero of the many-body 
wavefunction and binds a single flux quantum. These full 
QHs will be supplemented by a number of non-abelian 
fractional quasiholes at the appropriate filling fractions. 
In the table, each unitary braid matrix M is denoted by 
a pair of phases (pi,p2), where e ?,rpi and e mp2 are the 
eigenvalues of M. For cases with more than 2 impurities, 
we label the exchange of impurities i and j (as labeled in 
Fig. 1) by Rij. 

The simplest case v c r = 1/2 provides an excellent 
test of the algorithm, since we know (in the absence 
of a perturbing potential) that both the ground state 
wavefunction, and its quasihole excitations, are given ex- 
actly by Laughlin's variational ansatz I35| . On the torus 
the ground state is twofold degenerate [39j. Excitations 
about these two degenerate ground states states require 
overcoming an energy gap A ~ Jnn- The quasiholes are 
abelian anyons, and the Berry matrix in the ground state 
subspace should be the identity times a phase of ±tt/2, 
depending on the direction of the exchange path [jj 0] ■ 
This is consistent with our numerical studies of the path 
in Fig. 1(a). Since a complete braid of one quasihole 
around another is equivalent to two exchanges, we find 
a phase of ir for the path 1(c). As expected, when we 
introduce more impurities, we find that near is c g = 1/2 
all braids commute. 

A generic potential splits the two-fold degeneracy of 
the ground state by a small energy e. We attribute these 
splittings to interactions between the quasiparticlcs when 
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FIG. 1: (Color online) (a)-(c) Exchange paths used to braid 
quasiparticles on various lattices. In each path, the impu- 
rities (shaded red) are incrementally moved along the seg- 
ments (1,2,3...) until they return to their starting positions, 
exchanged. The dashed box represents the periodic lattice 
boundary, (d)-(i) show the initial configurations of the im- 
purities for the 3- and 4- impurity braids. 

they are moved close to one another. By optimizing the 
shapes of the potential at each time step, we can make 
e < 0.02A for all points in the v e s = 1/2 braid. While 
largely irrelevant for v c tf = 1/2, this optimization can 
be crucial for producing sensible results near v e g = 1 
or 3/2. If the trajectory is traversed in a time T such 
that h/A <C T <C h/e these splittings have no physical 
effect. We therefore neglect them when calculating M. 
We expect that the splittings can be further reduced by 
using larger systems. Detailed graphs of our optimized 
potentials are shown in the supplemental information for 
this paper [io| . 

The physics near z/ c h = 1 and 3/2 is richer. At z^ c ff = 1 
for ~Ui small, all particles are in the lowest Landau level 
and the ground state \1/g has a large overla p |l3j l with the 
Moore-Read (M-R) Pfaffian state * M r 0431lllil, a 
state with non-Abelian excitations. We typically per- 
form our calculations using hard-core interactions, where 
mixing with excited bands is significant and the over- 
lap is smaller: K^mrI^g)! < 0.3. Despite the small 
overlaps, one expects that the ground state with hard- 
core interactions is adiabatically connected to the M-R 
state and should share topological invariants such as ex- 
change statistics. The M-R state is the exact ground 
state of a Hamiltonian with repulsive three body inter- 
actions (f/2 = 0, XJ<& > 0). As we expand on below, we 
find excellent agreement between calculations using the 
two and three body interactions. 

The M-R state is gapped and has two types of funda- 
mental vortex excitations. In addition to the full QHs de- 
scribed earlier, the M-R state has half-quasiholc (HQH) 
excitations, which bind half a flux quantum, partially ex- 
clude particles from their location, and are non-abelian 
Ising anyons 0]. Wavefunctions of the M-R type with 
2n HQHs are 2™ _1 -fold degenerate [l9| in the limit that 
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i?34 : (0.69,-0.69) 



TABLE I: The results of our numerical braiding studies. Here, 
N is the total particle number, JV^ is the total number of 
flux quanta, and Mm P impurity sites have a repulsive poten- 
tial applied. "GFS" refers to whether the degenerate pair 
of eigenstates are the ground (G), first excited (F) or sec- 
ond excited (S) states. The braids are each characterized 
by a unitary matrix with eigenvalues e l7Tpi , e l7TP2 — >• (pi,P2). 
The exchange paths are shown in Fig. 1, with Rij denoting 
the exchange of impurities i and j. The algebras in the non- 
abelian cases approximate those described in the text [23l. l4ll| ; 
cases labeled as ambiguous contain non-commuting paths but 
the transformations associated with these paths depended on 
the details of the path and/or did not match the analytical 
predictions. Due to finite size splitting, not all paths were 
accessible on all lattices; only paths which led to a sensible 
braid and which were stable against small changes in the im- 
purity strength Vj are quoted here. Hard core interactions 
(U2 = 00) were used in all cases except 4x4*, where we also 
used (U2 =0,t/3 = 00). These two interactions gave nearly 
identical results. 

all the HQHs arc far apart. Given that we use strong 
impurity potentials (Vj > Jnn) we expect each repul- 
sive impurity will bind a full QH and a half quantum 
vortex. Exchanging two HQHs performs a 7r/2 rotation 
within the degenerate subspace, and the rotations pro- 
duced by exchanging different pairs of HQHs do not gen- 
erally commute. In particular, for four HQHs, it was 
shown [23I |47| that in the appropriate basis and ignoring 
Abelian phases, the braids can be written as 



i?12 = i?34 — e 



= i?24 — e 



(4) 



To estimate the overlap of the unitary transforma- 
tions which result from our braids with the predic- 
tions of the analytical theories of Bose quantum Hall 
states, we use the matrix overlap measure (M\ : M2) = 
|tr(MiM2)|/2. This quantity is insensitive to overall 
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phases, and we consider two unitary matrices to be equiv- 
alent if |tr(MiM 2 t )|/2 = 1. 

For the case of N = 4, = 8 and N[ mp = 4 on the 
6x4 lattice (where two impurities need never be near- 
est or next- nearest neighbors in a braid), our numerical 
results are in remarkable agreement with eq. La- 
beling the analytical predictions by R and the numerical 
matrices M, we have (Ri 2 ,Mi 2 ) = (#34,-^/34) = 0.99, 
(#24, M 2i ) = 0.98 and (R 13 ,M 13 ) = 0.97. 

When impurities are allowed to approach more closely, 
however, the numerical results diverge from the analyti- 
cal predictions, and in many cases, the exchange of two 
strong impurities produces a rotation by it. We conjec- 
ture that this represents the exchange of two pairs of 
HQHs, which either do not sit directly on the impurities 
but move with them as they are exchanged, or experi- 
ence tunneling events when impurities move too close to 
one another. For the case of a 4 x 4 lattice with N = 4, 



N, h 



8 and W 



imp 



4, we obtained identical results 



when considering the ordinary hard-core two-body or a 
hard core three-body interaction, where the M-R state 
is the exact ground state. For N — 7 and = 10 on 
the same lattice (Fig.QJ), we consistently obtained rota- 
tions by (0.5 ± 0.03) n, but the matrices which resulted 
were not straightforwardly related to the analytical pre- 
dictions in eq. (j4j) , and depended strongly on the path by 
which a pair of impurities were exchanged. These results 
show that the precise relationship of the non-abelian vor- 
ticies to the impurities is subtle 42-45|. Further, they 



reveal that the Berry matrices can be strongly modified 
for paths where impurities come close together. Surpris- 
ingly, the degeneracies are not necessarily broken by these 
close approaches. 

Finally, near v c g = 3/2 (fig. [Tg,h,j), we obtained a 
result consistent with the predictions for a Fibonacci 
anyon theory |4l|. the effective theory of the Read-Rezayi 
state at k = 3 [ll|. Previous numerical studies of con- 
tinuum bosons in the LLL [l3j have found strong evi- 
dence for this state, a particularly exciting result since Fi- 
bonacci anyons are capable of universal topological quan- 
tum computing. Comparing our numerically derived ma- 
trices at (N = 7, Nj, = 8 and N- lmp = 3) with the results 



of Hormozi et al. [4lJ, we obtained (i?i3,Mi3) = 0.99 
and (i?23, -^23) = 0.90. However, for the exchange of 
impurities 1 and 2, we found two sensible paths (a) and 
(b); in path (a) impurity 3 was allowed to move dur- 
ing the braid and in (b) it was not. We found that 
(Ri2,M 12 (a)) = 0.93, but (R 12 , M X2 (b)) = 0.69 and 
(M 12 (a),Mi 2 (b)) = 0.46. We conjecture that this dis- 
agreement was due to tunneling events when the impuri- 
ties were only next-nearest neighbors. For the 8x2 and 
9x2 lattices, we obtained rotations of nearly 37r/5 as 
predicted, but the resulting matrices had little overlap 
with those predicted from the Fibonacci anyon theory. 



SUMMARY AND CONCLUSIONS 

In summary, we have numerically studied a realistic 
model, eq. (JlJ, which has anyon excitations at filling 
fraction v e g = 1/2, and non-abelian anyons at u e g = 1 
and 3/2 analogous to those in the Moore- Read and Read- 
Rezayi states. These results suggest adiabatic continuity 
between the states of our lattice model with hard-core 
interactions and those found purely in the LLL (T3 |. to 
which our model reduces in the limit of weaker on-site 
interaction. We have also shown that surprisingly small 
lattices can reproduce infinite-system predictions, with- 
out resorting to trial wavefunctions. This robustness is 
likely related to the topologically protected nature of the 
states, and is encouraging for future experiments. 

The most intriguing implication of our result is in 
quantum computation. In recent years, a wealth of the- 
ory @, E3, 0, S3, S3 has shown that the M-R state of 
electrons at v = 5/2 could be used to construct topologi- 
cally protected quantum memory and quantum comput- 
ing operations, and has described potential implemen- 
tations. While non-abelian statistics in the v = 5/2 
state have not yet been confirmed experimentally, the 
fact that the v = 1 M-R state and the v = 5/2 M- 
R state are in the same universality class implies that 
the theory for manipulating quasiholes in the v = 5/2 
electron gas can be applied directly to our lattice bo- 
son system. Our v = 3/2 results are even more exciting 
since the Read-Rezayi states can be used to construct 
a universal quantum gate set. Implementing our model 
in a Josephson junction array would open a new area of 
physics to study topological noise protection and non- 
abelian statistics, since for <f> < 1/4 three non-abelian 
plateaux (y = 1,3/2, and 2) could be studied in the 
same experiment. The ability to individually address 
any lattice site would provide an unprecedented ability 
to manipulate quasiholes [IH , potentially creating a truly 
universal "quantum loom." 
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